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The morphological variation of basidiomata of Filoboletus manipularis (Berk) Singer collected in southern Viet Nam was 
studied. Phylogenetic analyses comprising three gene loci indicated that these collections, although exhibiting widely 
varying morphologies, represented a single species with a population composed of genetically diverse, sexually compatible 
monokaryon parental strains. No correlation was found between any aspect of morphological variation and intraspecific 
phylogenetic patterns for the three gene loci studied. Primers were designed to amplify the intron-rich 5' region of the 
translation elongation factor 1-a gene (tefla) and amplicons cloned and sequenced to characterize the parental haplotypes 
for individual basidiomata. The presence of recombination over the entire morphological diversity seen was confirmed by 
split decomposition analysis and analysis of gene diversity indicated a lack of allelic fixation within local populations. On 
several occasions, more than two apparent parental haplotypes were characterized from individual basidiomata, indicating 
that at least some basidiomata are chimeric or otherwise develop from a multinucleate condition. The literature supporting 
our observations of the occurrence of multinucleate basidiomata is reviewed and possible mechanisms for this phenomenon 
are proposed. 

Keywords: Filoboletus manipularis; morphology; population genetics; basidiomycete mating system; multinucleate 
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Introduction 

Filoboletus manipularis (Berk.) Singer (=Favolaschia 
manipularis (Berk.) Teng) is a luminescent fungal species 
widely distributed on rotting wood throughout Asian, 
Australasian and Pacific tropical regions (Corner 1954; 
Pegler 1986; Manimohan and Leelavathy 1989; Liu and 
Yang 1994; Desjardin et al. 2008). Fifteen species are 
currently accepted in the genus Filoboletus (http://www. 
mycobank.org); all described prior to 1996, with species 
circumscribed on the basis of morphological and anatomi- 
cal characteristics (e.g. Liu and Yang 1994; Maas 
Geesteranus and Horak 1995). The phylogenetic relation- 
ships among species accommodated in Filoboletus have 
not been studied, and even their generic affinities remain 
largely unresolved. Mycobank (http://www.mycobank.org/ 
) gives Favolaschia as the currently accepted generic 
placement for the species. However, phylogenetic analysis 
of available sequences (not shown) place F. manipularis 
within a Mycena clade closest to M. rubreomarginata and 
distant from Favolaschia species. According to Desjardins 
et al. (2008), F. manipularis and related species currently 
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assigned to Filoboletus and Poromycena require a new 
generic name since none are closely related to the type 
species of these genera. Publically available molecular 
data for Filoboletus is limited to five Genbank sequences 
that are fragments either of the 28S ribosomal DNA gene 
{Poromycena sp. AF261421, F. gracilis AF261422, 
Poromycena manipularis AF261423; Moncalvo et al. 
2002) or internal transcribed spacer 2 (F. manipularis 
AB509828, F. aff. manipularis AB509539). Therefore, 
the current study was undertaken to determine if the 
morphological variation observed in F. manipularis repre- 
sents more than one species or could be correlated to the 
infraspecific phylogenetics of the species. Collections of 
this species taken from southern Viet Nam in the current 
study exhibited a wide range of variation in morphologies, 
similar to morphological observations from previous stu- 
dies of the species from diverse tropical regions (Corner 
1954). DNA was isolated from individual basidiomata 
representing the diverse morphology observed and their 
genetic relationships determined from phylogenetic ana- 
lyses of three gene loci. 
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Methods 

Collections 

Individual and clustered basidiomata identified as F. 
manipularis were collected in Cat Tien National Park 
and Dong Nai Culture and Nature Reserve in southern 
Viet Nam (Table 1). These two reserves are core zones 
of the Dong Nai Biosphere Reserve and together com- 
prise 172,223 hectares of water surface (Tri An reser- 
voir), seasonally flooding grasslands and semi-evergreen 
and deciduous lowland tropical forests located in Dong 
Nai, Lam Dong and Binh Phudc provinces in Viet Nam. 
Collections were taken in June 2010, during the rainy 
season (average annual temperature and precipitation 
26°C and 2400 mm, respectively), in three separate 
regions of the reserves, near Nam Cat Tien, Bau Sau 
and Ma Da-VInh Ciru, located approximately 15-50 km 
distant from one another and representative of the habi- 
tat diversity of the forested portion of the reserves. All 
of the fruit bodies were found on the forest floor on 
rotting wood of undetermined origin, mostly in forests 
dominated either by Lagerstroemia calyculata Kurz or 
Dipterocarpus dyeri Pierre. Collected fruit bodies were 
air-dried on location and the separate collections 



transported for subsequent study in sealed plastic bags 
containing silica gel. Observations on the color (accord- 
ing to Kornerup and Wanscher 1978) and macromor- 
phology of basidiomata were recorded on fresh 
specimens in the field, and microscopic observations 
performed on the air-dried specimens as described pre- 
viously (Morozova et al. 2012) after transport to the 
investigating laboratory in Saint Petersburg, Russia. 
All specimens are deposited in the Komorov Botanical 
Institute Herbarium, Leningrad (LE). 

DNA isolation and sequencing 

Fruit body tissues were extracted in an aseptic manner 
from the pileus of individual basidiomata. Genomic 
DNA was isolated from the harvested tissues of each 
single basidiomata using the E.Z.N.A Fungal DNA Mini 
Kit (OMEGA Bio-Tek, Norcross, GA, USA) according to 
the manufacturer's protocol. The following genetic loci 
were amplified using the primers and amplification proto- 
cols listed in Table 2: (a) a region of nuclear rDNA 
containing the internal transcribed spacer regions 1 and 2 
(ITS1-2) and the 5.8S rRNA gene, (b) a fragment of the 
second largest subunit of the RNA polymerase II gene 



Table 1 . Collections of F. manipularis studied. 



Genbank no. 



Basidiome 


Collection no. a 


Date of collection 


Locality 


ITS 


rpb2 


tefla b 


Gl 


LE254356 


8 June 2010 


Nam Cat Tien 


KF746989 


KF746986 


KF746965 (6) 
KF746964 (5) 


G2 


LE 253916 


8 June 2010 


Nam Cat Tien 


KF746998 


KF746986 


KF746972 (9) 
KF746971 (4) 


G3 


LE 254357 


21 June 2010 


Ma Da-VInh Cuu 


KF746994 


KF746986 


KF746980 (7) 
KF746978 (2) 
KF746979 (2) 


G6 


LE 254358 


8 June 2010 


Nam Cat Tien 


KF746995 


KF746987 


KF746982 (9) 
KF746981 (2) 


G7 


LE 253912 


15 June 2010 


Nam Cat Tien 


KF746997 


KF746986 
KF746996 


KF746983 (19) 
KF746984 (4) 


G9 


LE 254359 


8 June 2010 


Nam Cat Tien 


KF746992 


KF746987 


KF746985 (9) 


Gil 


LE 254360 


12 June 2010 


Bau Sau 


KF746991 


KF746987 


KF746962 (2) 
KF746963 (2) 


G24 


LE 253914 


4 June 2010 


Nam Cat Tien 


KF746988 




KF746967 (6) 








Bau Sau 






KF746966 (3) 


G25 


LE 253911 


15 June 2010 


KF746990 


KF746987 


KF746969 (7) 
KF746970 (7) 
KF746968 (2) 


G30 


LE 253913 


15 June 2010 


Nam Cat Tien 


KF746993 


KF746986 


KF746975 (8) 
KF746974 (5) 
KF746973 (3) 


G31 


LE 253915 


4 June 2010 


Nam Cat Tien 






KF746977 (5) 
KF746976 (2) 



Sato et al.° 



Yakushima Island, Japan AB509828 



Notes: a LE - mycological herbarium of the Komarov Botanical Institute, RAS, St Petersburg. 
b in brackets - number of identical tefla clone sequences. 

c Sato H, Tsujino R, Kurita K, Yokoyama K, Agata K, 'Use of ITS sequences for resolving species richness, ecological differentiation and geographical 
distribution of macro-fungi in Yakushima Island, Japan', Genbank accession AB509828, 19 June 2009. 
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Table 2. 


PCR primers and protocols. 






Locus 


Primers 


Amplification protocol 


Reference 


ITS 

rpb2 
tefla 


ITS 1 -F (5 '-CTTGGTC ATTTAGAGGAAGTAA-3 ') 
ITS4-B (5'-CAGGAGACTTGTACACGGTCCAG-3') 

RPB2-filoF (5 '-ATTCCCCG ACC ATAACC AG-3 ') 
RPB2-filoR (5 '- AC AC AGGATGGCGACGATAG-3 ') 

TEF-filoF (5 '-TTTCTTGCAGCGCTTGTTCT-3 ') 
TEF-filoR (5 '-GTGCCAATACCACCGATCTT-3 ') 


3 min at 95°C; 35 cycles of (30 sec 

at 95°C, 30 sec at 55°C, 1 min at 

11 C); 10 min at 11 C 
3 min at 94°C; 40 cycles of (1 min 

at 94°C, 1 min at 65°C, 1 min at 

72°C); 10 min at 72°C 
30 sec at 98°C; 35 cycles of (10 sec 

at 98°C, 30 sec at 65°C, 30 sec at 

72°C); 10 min at 72°C 


Gardes and Bruns (1993) 
Current study 
Current study 



(rpb2) comprising approximately 250 bp and (c) a frag- 
ment at the 5' end of the translation elongation factor 1-a 
gene {tefla, eEFla) approximately 850 bp containing four 
introns. 

For ITS and rpb2, amplifications were performed 
using a Labnet Multigene thermocycler (Labnet Intl., 
Edison, NJ, USA). PCR reactions were prepared in 
10 uL volume containing the following mix: 1 uL 10X 
Titanium Taq Buffer (Clontech, Mountain View, CA, 
USA), 0.5 uL 2 mM dNTP, 0.32 uL each of 5 uM upper 
and lower primers, 7.76 uL sterile distilled water and 
0.1 uL Titanium Taq Polymerase (Clontech). For tefla, 
amplifications were performed using a Biometra 
TProfessional thermocycler (Biometra GmbH, 
Goettingen, Germany) with the Phusion Hot Start II 
High-Fidelity DNA Polymerase according to manufac- 
turer's protocol (Finnzymes, Espoo, Finland). 

Sequencing reactions were prepared using the ABI 
Prism® BigDye™ Terminator reaction kit (v3.1, Applied 
Biosystems Inc., Foster City, CA, USA) in 10 uL volume 
and 1/8 dilution using 5X sequencing buffer. The cycle 
sequencing reaction contained the following mix: 
1.75 uL 5X Sequencing Buffer, 0.5 uL BigDye V3.1 
Mix, 0.5 uL of 3.2 uM primer, 6.25 uL sterile distilled 
water, 1.0 uL (10-40 ng) PCR template, and employing 
the following amplification protocol: 25 cycles each of 
30-sec denaturation at 96°C, 15 sec annealing at 50°C 
and 4 min extension at 60°C. Primers used for PCR were 
also used for sequencing. Sequences were obtained using 
an ABI Prism 3100 Genetic Analyzer (Applied 
Biosystems). 

Cloning was performed on amplicons arising from 
amplification of tefla using the TOPO TA Cloning® Kit 
(Invitrogen, Carlsbad, CA, USA) following the manufac- 
turer's protocol. PCR for tefla was performed using a 
high-fidelity polymerase to reduce the frequency of clon- 
ing errors during nucleotide transcription. Nevertheless, 
12 singleton tefla clone sequences were excluded from 
phylogenetic analyses and only clone sequences obtained 
two or more times from individual basidiomata were 
included to reduce the possibility or frequency of cloning 



errors. Genbank accession numbers for the sequences are 
provided in Table 1. 



Phylogenetic analysis 

DNA sequences were aligned using the multiple sequence 
alignment program MAFFT and G-INS-i strategy (http:// 
mafft.cbrc.jp/alignment/server/), or the FFT-NS-i strategy 
for tefla fragments containing introns. The coding region 
and introns in the tefla locus were located from an alignment 
with annotated sequences of Mycena plumbea (Genbank 
GUI 87729) and Phyllotopsis sp. (DQ059047) and introns 
were removed for analyses of the tefla locus that included 
outgroup taxa. Parsimony phylogenetic analyses to visualize 
the degree of genetic variation for all three loci were per- 
formed with PAUP 4.0b. 10 (Sinauer Associates, Inc., 
Sunderland, MA, USA; Swofford 1993), using a heuristic 
search with a starting tree obtained via stepwise addition, 200 
Maxtrees and otherwise default parameters. Bayesian infer- 
ences were used for phylogenetic analyses of the tefla 
region. Sequences containing introns and with the intron 
regions deleted were analyzed separately, the latter analysis 
including Mycena plumbea (GUI 87729) as the phylogeneti- 
cally closest available outgroup representative determined 
from BLAST comparison of GenBank sequences. 
Nucleotide substitution models were determined using 
MrModeltest 2.3 (Nylander 2004; http://www.softpedia. 
com/get/Science-CAD/MrModeltest.shtml). The uncon- 
strained substitution model (GTR + I + G) was selected by 
both Akaike and Bayesian Information Criteria and applied 
in all analyses of tefla fragments. Bayesian estimations of 
phylogeny based on Metropolis-coupled Markov chain 
Monte Carlo sampling were performed using MrBayes 
v.3.2.1 (http://mrbayes.sourceforge.net/), with two simulta- 
neous runs of four incrementally heated chains running for 
ten million generations. Posterior probabilities were obtained 
from 50% majority rule consensus trees sampled every 5000 
generations. Probability values for posterior probabilities 
greater than 90% are shown on the phylograms. Patterns of 
recombination were visualized by split decomposition in 
SplitsTree v.4.12.3 (Huson and Bryant 2006; http://www. 
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splitstree.org), using pairwise distances with the Kimura 
K3ST model. The presence or absence of recombination 
was determined using the pairwise homoplasy index (PHI) 
implemented in SplitsTree. The level of heterozygosity 
among individuals within populations (localities) and 
among populations was determined from the tefla exon 
data using the program FSTAT to analyze gene diversity 
(http://www2.unil.ch/popgen/softwares/fstat.htm). 



Results 

Morphological observations 

Considerable variation was seen in the size, shape, color 
and luminosity of basidiomata (Table 3, Figures 1-6), 
although fruiting bodies arising within a single cluster 
were relatively uniform in morphology. At maturity fruit- 
ing bodies varied from 0.5 to 6.0 cm in pileus diameter 
with stipes 2.0 to 7.0 cm in length. The pileus could be 
conical, rounded, plane or depressed and umbonate or not. 
Color at maturity varied from white to cream or beige 
(near 4A2 'yellowish white', Kornerup and Wanscher 
1978) or pale pink (6A2 'orange white'). Developing 
basidiomata initially could be white or brownish (5C4 
'brownish orange') with the intensity of brown pigment 
diminishing with maturity. There were conspicuous differ- 
ences in luminosity patterns with only the cap luminescent 
or luminescent from underneath, only the stipe luminous, 
the entire fruiting body luminescent, or showing no lumi- 
nescence throughout. 

Microscopic observations indicated a degree of corre- 
lation between the overall gross morphology and the 
microscopic anatomy of the basidiomata. In most 
instances larger basidiomata had noticeably larger basidia, 
shorter and mostly simple-lageniform cheilocystidia, 
much larger caulocystidia and pileocystidia absent or 
nearly so. Smaller basidiomata had cheilocystidia mostly 
diverticulate or irregularly branched, less ornate caulocys- 
tidia, and abundant, irregularly shaped pileocystidia. 
These differences allowed the differentiation of the two 
reasonably distinct morphologies described below. 

Filoboletus manipularis (Berk.) Singer, Lloydia, 8: 215, 
1945 - Favolus manipularis Berk., Hooker's Journal of 
Botany and Kew Garden Miscellany, 6: 229, 1854; 
Laschia caespitosa var. manipularis (Berk.) Sacc, 
Sylloge Fungorum, 6: 407, 1888; Laschia manipularis 
(Berk.) Sacc, Sylloge Fungorum, 6: 408, 1888; 
Poromycena manipularis (Berk.) R. Heim, Revue de 
Mycologie, 10: 35, 1945 ; Mycena manipularis (Berk.) 
Metrod, Prodr. Fl. Mycol. de Madagascar, III, Les 
Mycenes, 87, 1949, f. 51 [non M. manipularis (Berk.) 
Sacc, Sylloge Fungorum, 5: 272, 1887] ; Favolaschia 
manipularis (Berk.) Teng, Cung-kuo Ti Chen-chun, 
[Fungi of China]: 760. 1963. 



Morphology type 1 - Figure 7 

Pileus 0.5-2.0 cm in diam., conico-campanulate to con- 
vex with conical umbo, hygrophanous, translucently reti- 
culate, finely white pruinose, whitish, grayish to hyaline 
when moist, purely white to yellowish dried up. 
Hymenophore tubular, adnate or adnate-emarginate with 
a slightly decurrent tooth. Tubes 1.5-4 mm long, 
arranged in radial rows, numbering 5-7 in a row, with 
angular-round pores 0.5-1 mm wide, white. Stipe 20- 
60 x 0.5-2.5 mm, cylindrical, thickened in the base, 
hollow, white to hyaline, completely white pruinose. 
Odor indistinct. 

Spores 6.2-7.5 * 4.4-5.1 urn, Q = 1.3-1.6, white, 
smooth, ellipsoid to broadly ellipsoid, amyloid. 
Basidia 16.3-18.3 x 6.2-7.5 um, clamped, narrowly 
clavate, with sterigmata up to 6 urn long. 
Cheilocystidia 51.5-73.1 x 6.8-12.9 um, forming a 
sterile lamellae edge, lageniform, subcylindrical or sub- 
clavate, mostly diverticulate or irregularly branched in 
the apex, or, rarely, simple lageniform with more or less 
developed neck. Pleurocystidia not seen. Hyphae of the 
pileipellis 5-15 um wide, clamped, with rare excres- 
cences and diverticulate terminal cells. Pileocystidia 
25.5-32.9 x 8.2-10.8 um, abundant, irregularly shaped, 
lageniform to subclavate, with excrescences. Hyphae of 
the cortical layer of the stipe 3-8 um wide, clamped, 
with rare excrescences and diverticulate terminal cells. 
Caulocystidia 24.4-38.1 x 6.8-7.4 um, abundant, irre- 
gularly shaped. 

Representatives. Viet Nam, Dong Nai Prov., Tan Phu 
Dist., Cat Tien National Park, on fallen log in lowland 
semideciduous tropical forest, 4 June 2010, col. O. 
Morozova, LE 253914; 15 June 2010, LE 253912. 



Morphology type 2 - Figure 8 

Pileus 1.5-6.0 cm in diam., convex, plano-convex to 
depressed, sometimes with a small papilla in the central 
depression, hygrophanous, translucently reticulate, 
smooth, slightly lubricious when moist, whitish with a 
pinkish tinge, with reddish brown spots in the old basi- 
diomata. Hymenophore tubular, adnate, adnate-emargi- 
nate with a decurrent tooth or shortly decurrent. Tubes 
2-5 mm long, adnate or adnate-emarginate with a 
slightly decurrent tooth, in radial rows, with angular 
pores 0.5-1.5 mm wide. Stipe 20-70 x 0.5-2.5 mm, 
cylindrical, thickened in the base, hollow, white with a 
pinkish tinge, brownish toward the base, completely 
white pruinose. 

Spores 5.8-8.6 x 4.0-5.5 um, Q = 1.3-1.9, white, smooth, 
ellipsoid to broadly ellipsoid, amyloid. Basidia 17.4— 
35.1 x 7.4-9.7 um, clamped, narrowly clavate, with 
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Figures 1-6. Filoboletus manipularis from southern Viet Nam showing variation in basidioma gross morphology. Figure 1: LE 254357 
(G3), white basidiomata of medium size and conical-umbonate pilei, photograph A.E. Kovalenko. Figure 2: LE 254357 (G3) showing 
stipes luminescent from below, photograph A.E. Kovalenko. Figure 3: LE 253914 (G24) white basidiomata with umbonate pileus, 
photograph O.V. Morozava. Figure 4: LE 253917 (G22) pink basidiomata with rounded pileus, photograph O.V. Morozava. Figure 5: LE 
253916 (G2) large, pink basidiomata with plane or depressed pileus. Figure 6 : LE 253912 (G7) cluster of cream to beige colored 
caespitose basidiomata. 



sterigmata up to 9.5 um long. Cheilocystidia 20.5- 
47.9 x 8.2-10.6 um, mostly simple lageniform or subcla- 
vate, rarely bifurcate in the apex. Pleurocystidia not seen. 
Hyphae of the pileipellis 5-15 um wide, clamped, with 
rare excrescences and diverticulate terminal cells, gelati- 
nized. Pileocystidia absent. Hyphae of the cortical 
layer of the stipe 3-8 um wide, clamped, with rare excres- 
cences and diverticulate terminal cells. 
Caulocystidia 20.5-138.0 x 6.5-11.0 um, abundant, irre- 
gularly shaped. 

Representatives. Viet Nam, Dong Nai Prov., Tan Phu 
Dist., Cat Tien National Park, on fallen log in lowland 
semideciduous tropical forest, 3 June 2010, coll. O. 



Morozova, LE 253917; 4 June 2010, coll. O. Morozova, 
LE 253915; 8 June 2010, coll. O. Morozova, LE 253916; 
21 June 2010, coll. A. Kovalenko, LE 254357. 

Other collections, however, showed more or less 
intermediate or mixed anatomical features not readily 
assigned to either morphology, and LE 253911 was 
notably divergent in having much smaller, narrow and 
irregularly shaped cheilocystidia, otherwise with basidia 
and pileocystidia similar to the type 1 morphology but 
with caulocystidia having the type 2 morphology 
(Figure 9). Observations that morphological variation 
was more or less continuous preclude the establishment 
of formal 'morphotypes' for the species based on these 
observations. 
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Figure 7. Filoboletus manipularis morphotype 1, LE 253914 
(G24) - characterized by small basidia, relatively large cheilo- 
cystidia and much reduced caulocystidia. (a) basidiospores, (b) 
basidium, (c) cheilocysidia, (d) pileocystidia and (e) caulocysti- 
dia. Bars = 10 urn. 



Phylogenetic analyses 

ITS was sequenced for 13 F. manipularis collections, 
with the ITS1-5.8S-ITS2 region varying from 667 to 
669 nt in length. Four sequence polymorphisms were 
found in the ITS regions from the sequence alignment 
(Table 4). These were an adenosine indel at position 
195 and a T-C transition at position 338 in the 1TS1 
region, and a second adenosine indel at position 538 
and A-G transition at position 669 in the 1TS2 region. 
The 5.8 region and bordering fragments of 28S and 
18S were invariant in our samples. Two sequence var- 
iants were noted for collections G3, G7, G24 and G25 
and unbalanced indels occurred in sequences for G3, 
G24 and G25, indicating ITS sequence differences 
including length variants between haplotypes from 
individual basidiomata in these collections. In total, 
there were 10 ITS sequence variants noted for all of 
the collections and median-joining analysis of the poly- 
morphism similarities indicated an interconnected, non- 



branching genetic relationship among all presumed 
dikaryon strains (Figure 10). The ITS2 sequence for 
F. manipularis from Japan (Genbank AB509828) was 
identical to sequences for eight of the Viet Nam col- 
lections indicating a close genetic relationship over a 
broader geographic region than sampled in the current 
study. 

rpb2 

We were unable with Filoboletus to successfully apply 
published protocols to amplify rpb2 in basidiomycetes 
(Morehouse et al. 2003; Matheny et al. 2007). Herein we 
designed primers to amplify a fragment of rpb2 for F. 
manipularis. These primers succeeded for all of our 
Filoboletus collections, although they did not work for 
Mycena. However, the length of amplified rpb2 fragments 
was less than 250 bp, and this gene fragment was not 
sufficiently informative to further elucidate the genetic 
structure of the sampled population of F. manipularis. 
Only two sequence polymorphisms were seen in this 
gene fragment and both were T-C transitions occurring 
at positions 31 and 57 of the alignment (not shown) 
discriminating two sequence types. Both rpb2 sequence 
variants were observed in G16 reflecting sequence differ- 
ences between the presumed parental monokaryons of this 
basidioma. 

tefla 

Herein we designed primers to sequence the 5' region 
of tefla, which we expected to show a relatively high 
degree of variation due to the presence of introns. 
These primers amplified a tefla fragment 843-858 bp 
in length containing 195 codons interspersed with four 
introns, the latter varying in length from 50 to 58 bp. 
After trimming of the 5' and 3' ends and deletion of the 
intron regions, all sequences including Phyllotopsis sp. 
(Genbank DQ059047) and Mycena plumbea 
(GUI 87729) were uniformly 586 nt in length. 
Polymorphisms in the tefla exon regions are shown 
in Table SI (Supplementary data) and the frequencies 
of various types of allelic polymorphisms in F. manip- 
ularis are summarized in Table 5. All 32 polymorph- 
isms observed in the exon regions were single 
nucleotide polymorphisms (SNPs). The majority of 
SNPs were nucleotide transitions (27, 84%) with C-T 
transitions predominant (22, 69%), and with only five 
nucleotide transversions and no indels. In the intron 
regions, there were a total of 97 polymorphisms, 
including 51 (53%) nucleotide transitions, 34 (35%) 
transversions and 12 indels of 1-4 bp length. C-T 
transition polymorphisms, as in the exons, predomi- 
nated in the intron regions with 44 occurrences 
(45%). Overall, the frequency of polymorphisms in 
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Figure 8. Filoboletus manipularis morphotype 2, LE 253916 (G2) - characterized by comparatively large basidia, shorter and simpler 
cheilocystidia and much reduced pileocystidia. (a) basidiospores, (b) basidium, (c) cheilocysidia, (d) pileocystidia and (e) caulocystidia. 
Bars = 10 um. 



the intron regions (97/207 bp = 46.8%) was approxi- 
mately nine times higher than in the exons (32/ 
586 pb = 5.4%). Three or more nucleotide types were 
observed at 12 positions in the introns and twice in the 
exons. The high number of polymorphisms in the 
introns increases the likelihood of multiple mutations 
including potential reversals occurring at the same 
locus; therefore, diversity seen in the exon regions 
should be phylogenetically more informative. 

The phylogenetic analysis based on tefla exon 
regions with Mycena plumbea (GUI 87729) as outgroup 
(Figure 11, http://purl.org/phylo/treebase/phylows/study/ 
TB2:S14862) indicated considerable genetic distance 
between Filoboletus and Mycena with limited phyloge- 
netic structure within F. manipularis. Analysis of gene 



diversities showed an unusually low level of hetero- 
zygosity among populations from the three different 
localities within the study area (Fst = -0.023, 
a = 0.060, supplementary Table SI) indicating an 
absence of allelic fixation or inbreeding within loca- 
lized populations. In the phylogenetic analysis fifteen 
clone sequence variants were grouped together in a 
polytomy collapsed at the 50% significance level. 
Among these, one clone variant (i.e. haplotype repre- 
senting one presumed parental monokaryon) from G2 
had identical exon sequence to one haplotype from 
each of G7 and G30, and one haplotype from G3 
was identical to two haplotypes from G31. Similarly, 
Gl and G7 had one haplotype with identical exon 
sequences. However, none of the basidiomata had 
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Figure 9. Filoboletus manipularis LE 25391 1 (G25) - characterized by reduced or shrunken cheilocystidia. (a) basidiospores, (b) basidium, 
(c) cheilocysidia, (d) pileocystidia, (e) caulocystidia. Bars = 10 urn. 



identical paired haplotypes, indicating a rich assem- 
blage of compatible mating strains in relatively close 
proximity to one another in these locations. Only two 
branches within Filoboletus were supported at levels 
higher than 90% and none higher than 95%. One 
branch joined two haplotypes from G3, and the second 
joined haplotypes from Gl and G7. However, all three 
basidiomata had one additional haplotype outside of 
the supported branches. Significantly, three different 
haplotype sequences were seen for basidiomata of G3, 
G25 and G30 (supplementary Table SI). Since DNA 
for this study was isolated from individual basidiomata, 
which are presumed to have arisen from conjugation of 
two compatible monokaryons, this indicates that more 
than two haploid components could be involved in the 
development of individual basidiomata, 



Similar observations were found when the intron regions 
were included in the analysis (Figure 12, http://purl.org/ 
phylo/treebase/phylows/study/TB2:S 14862). Significant 
branches (>0.95) joined paired haplotypes for G31, G25 
and G3, as well as for individual haplotypes from Gl with 
G7, and G2 with G24. With the exception of G3 1, the other 
basidiomata had additional haplotypes arising on unsup- 
ported branches (<0.95), again indicating the possibility 
that more than two component nuclei were involved in the 
development of individual basidiomata. G7 and G30 each 
had one haplotype with identical sequence, as did Gl and 
G7; however, all of these collections originated from the 
vicinity of Nam Cat Tien. Recombination patterns were 
visualized by split decomposition analysis (Figure 13) onto 
which the morphological characteristics of the basidiomata 
were summarized. This analysis showed a complex network 
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Table 4. ITS polymorphisms. 



Polymorphism 



Isolates 


1 Q^ a 
LyJ 


1 1. 

JJO 


^1S C 


oov 


G2, G16, G19 


A 


T 


A 


G 


G25A 


A 


T 


A 


A 


G6, Gil, G13 


A 


T 




G 


G25B 


A 


T 




A 


Gl, G7A, G9, G24B 




C 




G 


G7B 




C 




A 


G24A 




c 


A 


G 


G3A 




T 


A 


G 


G30 




T 




A 


G3B 




T 




G 


AB509828 








G 



Notes: "Adenosine indel at position 195 in ITS1. 
T-C transition at position 338 in ITS1. 
"Adenosine indel at position 538 in ITS2. 
d A-G transition at position 669 in ITS2. 
- Missing sequence for ITS 1 region. 
■ nt deletion. 



G24A 



G3A 



G1, G7A 
G9, G24B 




G7B 



G30 



Figure 10. Median-joining network, showing polymorphisms 
in the internal transcribed spacer (ITS) regions. Lines separat- 
ing ITS haplotypes represent polymorphisms: (a) adenosine 
indel at position 195, (b) T/C transition at position 338, (c) 
adenosine indel at position 538, (d) A/G transition at position 
669; * alternate possible positions for AB509828 from Japan. 



of relationships and statistically significant recombination 
within all of F. manipularis (Phitest, <t> = 3.64E-11). There 
was no evident pattern of phylogenetic differentiation or 
correlation of genetic variation with morphological 
characteristics. 



Table 5. Number and relative frequency of allelic polymorph- 
isms in the tefla regions of F. manipularis. 



Mutation 


Exons 


Introns 


Total 


A-G 


5 (16%) 


7 (7%) 


12 (9%) 


C-T 


22 (69%) 


44 (45%) 


66 (51%) 


Total transitions 


27 (84%) 


51 (53%) 


78 (60%) 


A-C 


2 


9 


11 


A-T 


1 


12 


13 


C-G 


2 


6 


8 


G-T 


0 


7 


7 


Total transversions 


5 (16%) 


34 (35%) 


39 (30%) 


Indels 


0 


12 (12%) 


12 (9%) 


Total mutations 


32 


97 


129 


Length (nt) 


586 


207 


793 



0.1 



— G11A 2 clones 
G1A 5 clones 
-G31A 2 clones 

- G31B 5 clones 

- G3A 7 clones 

- G2A 9 clones 

- G30A 3 clones 

- G7A 19 clones 
G24B 6 clones 
G6A 9 clones 
• G30B 5 clones 

- G2B 4 clones 
G25A 7 clones 

G9 9 clones 

G6B 2 clones 

G1B 6 clones 
G7B 4 clones 
G3B 2 clones 
G3C 2 clones 

- G11B 2 clones 
G24A 3 clones 

— G30C 8 clones 

G25B 7 clones 

G25C 2 clones 




GU1 87729 Mycena plumbea 



Figure 1 1 . Fifty percent consensus tree from Bayesian analysis 
of tefla coding region (exons). Posterior probabilities >0.90 
indicated above statistically significant branches. 



Discussion 

No apparent phylogenetic basis for observed 
morphological variation in F. manipularis 

Considerable variation was observed in basidioma mor- 
phology among the Viet Nam collections, similar to the 
variation previously reported on a worldwide scale 



(Corner 1954; Pegler 1986; Manimohan and 
Leelavathy 1989). Basidiomata at apparent maturity 
varied in size (0.5-6.0 cm diam.), shape (conical, cam- 
panulate or plano-convex to depressed) and color (white 
to cream, beige or pinkish). Similarly, patterns of lumi- 
nescence were highly variable, with either the pileus or 
stipe luminescent, or the entire basidioma luminescent 



Mycology 9 1 



G25A 7 clones 

G30C 8 clones 




G3C 2 clones 
G3B 2 clones 

G3A 7 clones 

G6B 2 clones 

G6A 9 clones 
G24A 3 clones 



G3 1 A 2 clones 
G31 B 5 clones 

G11B 2 clones 

G25C 2 clones 
7 clones 

G30B 5 clones 



Figure 12. Bayesian analysis of tefla locus including introns. 
Open circles indicate posterior possibilities (PP) > 0.90; closed 
circles PP > 0.95; branches collapsed where PP < 0.50. 



or not. Notably, larger basidiomata were often pink- 
tinged and appeared to have larger basidia, smaller, 
essentially inornate cheilocystidia and much larger cau- 
locystidia than smaller basidiomata which tended 
toward white at maturity. The collection corresponding 
to G25 (LE 253911) was aberrant in having abnormally 



small, narrow cheilocystidia with conspicuous lateral 
invaginations as though collapsing. These are anatomi- 
cal distinctions that have frequently been used to dis- 
tinguish species based on morphology. However, we did 
not find any correlates between any aspect of morpho- 
logical variation and intraspecific phylogenetic patterns 
for any of the three gene regions studied. Luminosity 
also was uncorrelated to any aspect of morphological 
variation or phylogenetic pattern. Overall, the variation 
described by Corner (1954) for collections from diverse 
tropical regions essentially matches our observations. 
Corner had the opportunity to observe the development 
over time of fruit bodies of F. manipularis collected in 
northeastern Pakistan. He noted that the fruit bodies 
were slow in their expansion and remarkably long- 
lived, producing spores variably for 2-6 days. 
Significantly, he was able to observe age differences in 
the appearance of cheilocystidia, caulocystidia and 
pileocystidia, with cheilocystidia and pileocystidia col- 
lapsing on expansion of the pileus. This provides a 
likely explanation for the aberrant and apparently senes- 
cent cheilocystidia observed in LE 253911 in our study. 
Corner (1954) also noted a progression of color changes 
over the life of the basidiomata, and concluded that pink 
and reddish discoloration in age might be attributable to 
the presence of molds (e.g. Poldmaa 2011). 
Conceivably, genetic variation between sexually compa- 
tible monokaryons could contribute to the 



Basidiome color 
W white 

P pink Luminescence 

□ none 
H pileus 

□ stipe 
■ all 




' G25A 7 clones 
IG7A 19 clonesV r\ 

/G30A 3 clonesS n 

/ Q Q G24B 6 clonesH , 

■ □ o_ 



s 



w 



w 



G1B 6 clones G3C 2 clones 
Q G3B 2 clones 

□ 
W 



□ 



Basidiome size 

O small < 25mm diam 

% large > 25mm 



Figure 13. Split decomposition analysis of the tefla locus. Pairwise homoplasy index (PHI) = 3.646 , indicating the presence of 
recombination across all collections. 
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morphological variation observed. However, environ- 
mental factors such as climate, local weather patterns 
or substrate suitability, and developmental factors such 
as age in particular are very likely to be significant 
determinants of basidioma morphology and luminosity. 
In the current study anatomical observations may have 
been significantly impacted by varying conditions of 
drying or preservation since the collections had to be 
shipped and were not immediately available for micro- 
scopic observation. All these factors advise caution in 
interpreting morphological and anatomical characters to 
differentiate and define species in Filoboletus. 

Filoboletus manipularis is a single phylogenetic species 
as currently recognized 

Only four polymorphisms were found in a 667-669 bp 
sequence fragment of ribosomal DNA which included the 
ITS 1 and 2 regions. This represented a level of ITS 
variability of 0.6%, which is well within the generally 
accepted level of intraspecific ITS variation of approxi- 
mately 0-3% (Basidiomycota average 3.3%, Nilsson et al. 
2008). The only published sequence for F. manipularis 
(Genbank AB502898 from southern Japan) also fell 
within this range of variation having an identical ITS2 
sequence to 6 of the 10 Vietnamese collections. Unique 
haplotype sequences were observed in G7 and G30. All 
other collections shared identical haplotype sequences in 
varying combinations, indicating a complexly intercon- 
nected relationship among all of the collections. In con- 
trast, 'Filoboletus aff. manipularis'' (AB509539), when 
compared to F. manipularis (AB502898), had 24 poly- 
morphisms within the 1TS2 region (209 nt's) including 9 
indels of 1-2 bp, thus representing a frequency of ITS 
variation of 11.5% and indicative of a separate species. 
Our observations of an approximately 250 bp fragment of 
rpb2 also indicated a close, interconnected phylogenetic 
relationship for the collections studied, with only three T- 
C nucleotide transitions differentiating the collections. 
Identical haplotype sequences also were observed from 
morphologically distinct basidiomata in analyses of the 
tefla locus, reinforcing the conclusion that the Viet Nam 
collections represent a single species. The estimated coef- 
ficient of inbreeding was unusually low (Fst = -0.023), 
indicating a lack of allelic fixation within local popula- 
tions sampled in this study. Conclusively, the analysis of 
split frequencies based on the tefla sequence data found 
statistically significant recombination among all of the F. 
manipularis collections and no correlation between 
genetic relationships and morphological characteristics. 

Could basidiomata be 'polykaryotic' or chimeric? 

The detection of three different tefla haplotype sequences 
from single basidiomata of G3, G25 and G30 indicates a 



possibility that more than two monokaryons contributed to 
the formation of individual basidiomata. A high fidelity 
'proofreading' polymerase was used for PCR after clon- 
ing; nevertheless, 12 unique ('singleton') clone sequences 
were excluded from the phylogenetic analyses to eliminate 
the possibility of cloning errors caused by the incorrect 
incorporation of nucleotides. Assuming the occurrence of 
mutations is random, the possibility of any two sequences 
from one basidioma having the same mutation can be 
approximated using the calculation provided by 
Vydryakova et al. (2011) based on the 'Birthday 
Paradox'. For G3, G25 and G30, there are 15-17 poly- 
morphic regions and the possibility of random variation 
explaining two occurrences of any one transition poly- 
morphism is approximately P ~ 0.148 and of three occur- 
rences P ~ 0.022, or half these probabilities for a 
transversion mutation. For G3 the two 2-clone sequences 
differ by a single base, a G-C transversion, which could 
statistically (P ~ 0.074) be attributed to cloning error, that 
if correct would leave two haplotypes represented by 7 
and 4 clones, respectively, and differing by 17 distinct 
polymorphisms. For G25, there are 28 total polymorph- 
isms differentiating the 2-clone sequence and the two 
7-clone sequences. However, the 2-clone sequence has 
only two unique polymorphisms, otherwise having one 
or the other polymorphisms of either of the two 7-clone 
sequences, suggesting it could be a recombinant derived 
from the 7-clone sequences with one crossover occurring 
between positions 513-579 of the alignment (Table 6). 
However, the probability of occurrence of the two unique 
polymorphisms in the 2-clone sequence, a C-T transition 
and an A-T transversion, is P-0.011. Therefore, it would 
appear more likely that the 2-clone sequence is a unique 
third haplotype in G25. For G30 there are 29 polymorph- 
isms involving the 3-, 5- and 8-clone sequences (Table 7). 
The 3 -clone sequence in G30 is unique in four polymorph- 
isms, including a 4-base substitution, two C-T transitions 
and a 2-base deletion that statistically cannot be attributed 
to cloning error leaving the additional possibility of three 
parental monokaryons contributing to the G30 basidioma. 

We had a closer look at the 12 unique clone sequences 
('singleton' clones) excluded from the previous analyses 
to see if they showed significant patterns of variation. Four 
clones had unique polymorphisms that could have resulted 
from PCR replication errors propagated by cloning. All 
polymorphic loci in the remaining eight singleton clone 
sequences had the same polymorphisms as one or more 
other clones. Many clones had polymorphisms in common 
with one or another of the dominant clones suggestive of 
recent recombination involving the dominant clones. In 
G3, the pattern of polymorphisms in all four singleton 
clones could be explained by one or two crossover events 
occurring within a span of about 700 nts (Table 8). 
However, the singleton clone in G6 (Table 9) would 
require eight crossover events to have been derived by 
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Table 6. tefla polymorphisms in G25 clones. 



G25 clones 


34 


40 


156 179 194 207 219 


222 


Locus a ' b 

341 419 456 461 466 493 496 510 513 551 579 649 657 660 668 672 676 738 765 768 


7 clones B 


GCCA 


C 


C ■ C C G 


TGCT 


C 


C 


CTATCCTACACC ■ ATCCA 


2 clones C 




T 


G C T A C 


GCTA 


T 


T 


■ CGCTTCT CACC ■ ATCCA 


7 clones A 




T 


G C T A C 


GCTA 


T 


C 


■ CGCTTCATC ■ TTCTCTTG 



Notes: a Intron loci italicized, exon loci in bold type, start position for consecutive polymorphisms, indels. 
b Shared polymorphisms shaded, possible crossing-over location indicated by vertical lines. 
• nt deletions. 



Table 7. tefla polymorphisms in G30 clones. 



G30 clones 


67 


168 


ISO 


184 


194 


208 


219 


222 


245 


416 


440 


453 


456 


Locus* 
460 466 


472 


481 


485 


493 


579 


64< 


650 


656 


660 669 


744 


765 


768 


817 


5 clones B 


T 


C 


C 


T 


C 


C 


C 


GCTA 


C 


C 


T 


G 


11 


T G- 


CAGC 


a 


A 


T 


C 


A 


A 


C 


C 


T 


C 


G 


T 


3 clones A 


T 


C 


T 


T 


C 


C 


C 


GCTA 


C 


C 


C 


G 




T G- 


ATGT 


G 


A 


C 


C 


1 


C 


T 


T- 


C 


T 


G 


T 


8 clones C 


C 


T 


T 


C 


T 


G 


G 




T 


T 


C 


C 




C AC 


AACC 


T 


G 


C 


T 


C 


C 


T 


• TTT 


C 


C 


A 


C 



Notes: a Intron loci italicized, exon loci in bold type, start position for consecutive polymorphisms, indels. 
b Shared polymorphisms indicated by shading. 
• nt deletions. 



Table 8. tefla polymorphisms in G3 clones. 



G3 clones Locus ab 

49 51 67 179 245 350 455 459 463 474 510 643 648 655 657 659 691 753 817 



7 clones A 


T 


T 


C 


T 


T 


c 


C 


ATA 


T 


c 


C 


A 


C 


C 


T 


T 


C 


A 


c 


1 clone 


T 


T 


C 


T 


T 


c 


C 


ATA 


T 


c 


C 


A 


c 


c 


T 


T 


C 


C 


G 


1 clone 


C 


A 


T 


T 


T 


c 


c 


ATA 


T 


c 


C 


A 


c 


c 


T 


T 


C 


A 


C 


1 clone 


C 


A 


T 


T 


T 


c 


c 


ATA 


T 


c 


C 


A 


c 


c 


T 


T 


C 


A 


G 


1 clone 


T 


T 


C 




C 


T 


T 


TCT 


C 


G 


T 


A 


c 


c 


T 


T 


C 


A 


C 


2 clones B 


C 


A 


T 




c 


T 


T 


TCT 


c 


G 


T 




T 


T 


C 


C 


T 


C 


c 


2 clones C 


C 


A 


T 




c 


T 


T 


TCT 


c 


G 


T 




T 


T 


c 


C 


T 


c 


G 



Notes: a Intron loci italicized, exon loci in bold type, start position for consecutive polymorphisms, 
indels. 

b Shared polymorphisms indicated by shading, locations of possible crossings-over indicated by 
vertical lines. 
• nt deletions. 



Table 9. tefla polymorphisms in G6 clones. 



G6 clones 












Locus ab 




36 


51 


119 


207 


219 222 229 453 457 461 513 537 639 645 648 650 655 660 676 692 744 


9 clones A 


C 


A 


T 


ACT 


C GCTA A C • 


CCAACCCTC - CTC 


1 clone 


c 


A 


T 


ACT 


G • G G T 


C C G • T C C CA • T C C 


2 clones B 


T 


T 


C 


TTA 


G ■ G G T 


TTG ■ TTACACTCA 



Notes: a Intron loci italicized, exon loci in bold type, start position for consecutive polymorphisms, indels. 
b Shared polymorphisms indicated by shading, locations of possible crossings-over indicated by vertical lines. 
• nt deletions. 
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Table 10. tefla polymorphisms in G7 clones. 



G7 clones Locus a,b 

34 51 59 6 7 206 207 219 222 233 245 446 459 461 472 474 479 481 485 487 644 648 652 658 659 666 672 674 765 



19 clones A 


GCCA 


A 


T 


T 


C 


A 


c 


GCTA 


T 


C 


T 




T 


C 


G 


T 


G 


A 


C 


T 




C 


T 




G 


T 


T 


T 


1 clone 


GCCA 


A 


C 


C 


c 


A 


c 


GCTA 




C 


T 




T 


C 


G 


T 


G 


A 


C 


T 




c 


T 




G 


T 


T 


T 


1 clone 


GCCA 


A 


T 


C 


c 


A 


c 


GCTA 


T 


c 


T 




T 


C 


T 


T 


T 


G 


C 


T 




c 


T 


c 


G 


T 


T 


T 


1 clone 


GCCA 


A 


C 


C 


c 


A 


c 


GCTA 


C 


c 


T 




T 


c 


T 


T 


T 


G 


C 


T 




c 


T 


c 


G 


T 


T 


T 


1 clone 


GCCA 


A 


T 


C 


T 


A 


G 




T 


T 


T 


c 


c 


A 


G 


T 


T 


G 


C 


T 


T 


T 


C 




G 


T 


T 


T 


1 clone 


GCCA 


T 


T 


C 


c 


t 


G 




T 


T 


T 


t 


T 


A 


G 


T 


T 


G 


C 


c 


T 


c 


c 




G 




T 


T 


4 clones B 




T 


T 


T 


T 


G 


c 


GCTA 


T 


C 


C 


G 


A 


A 


G 


A 


T 


A 


T 


T 


C 


T 


c 




T 




A 


c 



Notes: a Intron loci italicized, exon loci in bold case, start position for consecutive polymorphisms, indels. 

b Pattern of shared polymorphisms indicated by shading, polymorphisms not shared with either predominant clone boxed, unique 
polymorphisms in bold and lower case. 
• nt indels. 



recombination from the dominant clones, with which it 
shares 2 1 polymorphic values in common with one or the 
other dominant clone. Clearly this level of crossing-over 
reflects a long evolutionary history of recombination 
events. A more complicated pattern of relationships is 
evident in the clones from G7 (Table 10). Singleton clones 
D and E have two and five unique polymorphisms, respec- 
tively, and three of these polymorphic loci are also poly- 
morphic between the dominant clones. In addition, 
polymorphisms occurring in differing combinations of 
singleton clones and not appearing in the dominant clones 
occur at ten polymorphic loci. Additional singleton clones 
were observed from basidiomes G30 and G31 (not shown) 
that were similarly indicative of haplotypes different from 
the dominant haplotypes. Overall, the pattern of distribu- 
tion of shared polymorphisms suggests the occurrence in 
individual basidiomata of multiple haploid components 
that have arisen following a long phylogenetic history of 
mutation and recombination. 

In planning this study, we believed we would find, at 
most, two different haploid sequences representing the two 
parental monokaryons contributing to formation of indivi- 
dual basidioma. However, more than two haploid sequences 
were found among clones from some basidiomata which 
could not be attributed to cloning errors arising from nucleo- 
tide misincorporation. Highly variable products of recombi- 
nation could occur in the generative tissues which could be 
propagated by PCR, and these recombinants could be com- 
mon if a significant portion of the genetic material is translo- 
cated from structural and generative tissues and devoted to 
sexual recombination. Patterns suggestive of recombination 
were seen for some of the 'singleton' clone sequences, as 
well as for two identical clones sequenced from G25. 
However, in most cases the pattern of recombination was 
overly complex to be attributable to recent recombination 
within a single basidioma, requiring frequent and close- 
spaced crossing-over. Finally, in G30 there were three dif- 
ferent and statistically significant clone sequences suggestive 



of three parental monokaryons, and multiple haploid 
sequences were also observed in G7 and G31 that could 
not be explained by recent meiotic recombination. 

Vydryakova et al. (2011) reported similar observations 
for the ITS locus in Neonothopanus nambi (Speg.) R.H. 
Petersen & Krisai, also suggestive of a 'polykaryotic' con- 
dition in basidiomata of this species. However, since ITS is a 
multicopy genetic locus these authors recognized the possi- 
bility that variability in the ITS could potentially be attrib- 
uted to the existence of multiple copies of the rRNA cistron 
evolving independently by being separately located in 
unlinked regions of the genome. Tefla has generally been 
considered to be a single-copy gene. However, there is 
increasing evidence for the occurrence of paralogous copies 
of the tefla gene in fungi (Aguileta et al. 2008). A search for 
the tefla gene in complete genomes of Agaricomycotina in 
Mycocosm using search term KOG0052 indicated that tefla 
is a single-copy gene in approximately 90% of the agaric 
species covered (http://genome.jgi.doe.gov/programs/fungi/ 
index.jsf, Grigoriev et al. 2012). Multiple tefla copies 
occurred in only 9 of the 94 genomes annotated to date, 
none of which were closely related to Filoboletus, and multi- 
ple copies were not seen in Tricholoma matsutake, which 
was the only representative of the Tricholomatoid clade 
(Matheny et al. 2006). The lack of full-genome coverage 
for the Mycenaceae leaves uncertainty concerning the prob- 
able copy state for the gene in Filoboletus. However, if 
multiple copies of tefla were amplified in the current study 
we would expect to see patterns of sequence divergence for 
the paralogous copies. Instead, the variation appears more or 
less random and continuous and split decomposition analysis 
of recombination patterns indicated a statistically significant 
level of recombination across all haplotype sequences. 
Therefore, we conclude that a single locus for the tefla 
gene was amplified in this study. An additional concern for 
the current study is the possibility of PCR-mediated recom- 
bination contributing to the sequence variation observed for 
tefla. Chimeras are produced when a prematurely terminated 
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amplicon reanneals to foreign DNA and is copied to com- 
pletion in a subsequent amplification. The high-fidelity poly- 
merase employed in this study requires a short extension 
period which has been associated with PCR-mediated 
recombination. However, this polymerase utilizes a low 
starting template concentration which should compensate 
to decrease the likelihood of chimeras being produced 
(Lahr and Katz 2009). It is difficult to distinguish chimeric 
sequences from sequence variants resulting from meiotic 
recombination. However, most of the haplotype sequences 
in the current study are distinguished by SNPs and indels in 
addition to crossovers that could have resulted from PCR- 
mediated recombination or meiotic recombination. Break 
points occurred at the exactly the same locations in multiple 
haplotypes and most haplotypes had a pattern of frequent and 
close-spaced crossovers unlikely to be seen in PCR- 
mediated chimeras. Instead the pattern of variation observed 
in this study is indicative of a long history of recombination 
and inherited mutations and provides compelling evidence 
for the occurrence of three or more parental monokaryons 
contributing to basidioma development. Two mechanisms 
that might result in multinucleate basidiomata are discussed 
below. 

Multinucleate basidiomata develop directly from 
trikaryotic or polykaryotic mycelium in multinucleate 
basidiomycetes 

Germinating basidiospores and mitotically generated oidia 
readily fuse with compatible mycelia to form the dikaryon 
mycelial phase that dominates the basidiomycete life 
cycle. However, the nuclear condition of the dikaryon 
phase is not stable, and di-mon mating - dikaryotization 
of a monokaryon following mating with an existing dikar- 
yon originally described by Buller (1930) - may occur 
frequently in nature. Stable trikaryons, following di-mon 
matings, have been observed in pairings between hetero- 
karyons and homokaryons in basidiomycetes with multi- 
nucleate cells (James et al. 2009). As many as 30% or 
more of basidiomycetes may have multinucleate cells 
('cenocyty', Boidin 1971). Therefore, if tri-genomic 
mycelia undergo recombination through completion of 
the sexual cycle (Kimura and Kadoya 1962), multinucle- 
ate basidiomata may occur relatively frequently in basi- 
diomycetes with multinucleate mycelium. Koltin and 
Raper (1968) also found evidence for the occurrence of 
a diploid mycelium in Schizophyllum capable of anasto- 
mosing with haploid mycelium to establish 'haploid- 
diploid dikaryons' in which the diploid nucleus appeared 
unstable resulting in a mosaic of two different 'haploid- 
haploid dikaryons' in the resulting basidiomata. However, 
in our observations clamp connections were formed reg- 
ularly in F. manipularis perhaps indicating against the 
possibility of a polykaryotic state being maintained to 
formation of basidiomata. 



Adjacent mycelium representing two different dikaryons 
conjoin to form a mycelial knot in the initial stage of 
basidiome formation with resultant chimeric basidiomata 

Following dikaryotization, hyphae readily fuse to allow 
active translocation of limiting nutrients throughout a 
mycelium (Olsson and Gray 1998). Nuclear replacement 
following anastomosis is unidirectional toward the hyphal 
tips (Kues 2000) and following sequential anastomoses 
originally separate clones may act as a single individual. 
Thus a mycelium may be a chimera of genetically differ- 
ent mycelial parts benefiting mutually in exploiting distant 
food sources by translocation of nutrients throughout the 
mycelial network. In addition, somatic recombination 
occurs at low frequencies in both trikaryon and dikaryon 
mycelium (Shalev et al. 1972; Korhonen and Hintikka 
1974; Frankel 1979; Anderson and Kohn 2007) and multi- 
ple nuclear haplotypes reassorted from somatically incom- 
patible mycelia have been found residing in the same 
mycelium (Johannesson and Stenlid 2004). Somatic 
incompatibilities do not necessarily result in the effective 
segregation of new-formed dikaryon pairings and a strictly 
compatible dikaryon condition is not always required for 
fruit body formation in basidiomycetes. Incompatibility at 
the B mating-type locus, for example, does not preclude 
fruit body formation and basidiomata can develop readily 
from common B heterokaryons (Casselton and Kues 1994; 
Kues 2000). Fruit body formation in the basidiomycetes is 
not necessarily monocentric and hyphal knots representing 
the initial stage in fruit body formation commonly origi- 
nate from more than one generative hypha, with branches 
from neighboring hyphae anastomosing to form the basi- 
dioma initial (Kues 2000, observations in Coprinus 
cinereus, = Coprinopsis cinerea, Psathyrellaceae). Since 
only the young hyphae at the margin of the colonies can 
produce basidiomata (Kues 2000), it is conceivable that 
adjacent mycelia with different nuclear types could come 
in contact at the time of fruiting. Genetic differences 
potentiate rejection and acceptance mechanisms; however, 
interaction outcomes are dependent on both genes and 
context and the ability of nonself genomes to conjoin 
may be dependent on levels of oxidative stress (Rayner 
1996). Similarly the frequency of cenocytic mycelium and 
multinucleate cells is increased under conditions of insuf- 
ficient aeration (Boidin 1971). The formation and toler- 
ance of chimeric mycelial communities through hyphal 
anastomoses advantage the supply of nutrition during 
fruit body development. In the absence of a strict nonself 
recognition mechanism, anastomoses between adjacent, 
genetically different dikaryons in the mycelial knot stage 
of basidioma formation could result in chimeric basidio- 
mata. Seamless chimeric basidiomata between genera 
have been recorded in two different situations. Kemp 
(1977) reported chimeric basidiomata jointly formed by 
Parasola misera (P. Karst.) Redhead, Vilgalys & Hopple 
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and Coprinellus pellucidus (P. Karst.) Redhead, Vilgalys 
& Moncalvo (both in Psathyrellaceae and reported as 
'Coprinus' species) that had colonized the same dung 
sample. The mycelia of the two species were clearly not 
antagonistic and he was able to distinguish the mycelium 
of the two species in the stipe as well as tetrads of spores 
of both species in the gills. Other chimeric basidiomata 
involve the mycoparasitic genus Squamanita. Mushroom 
species such as Squamanita contortipes (A.H. Sm. & D.E. 
Stuntz) Heinem. & Thoen and Squamanita paradoxa (A. 
H. Sm. & Singer) Bas, form chimeric basidiomata with 
their agaric hosts, a Galerina for S. contortipes, and a 
Cystoderma with S. paradoxa (Redhead et al. 1994; 
Matheny and Griffith 2010), as do other Squamanita spe- 
cies with their hosts (Reid 1983; Harmaja 1988; Nagasawa 
et al. 1990). 

In the case envisaged for Filoboletus, the apparent 
nonselfish collaboration of genetically different nuclei in 
the translocation of nutrients may serve to enhance the 
vigor of a chimeric mycelium allowing adaptation to 
changing or unfavorable environmental conditions. The 
existence of polykaryotic basidiomata also would be sig- 
nificant in reducing barriers to genetic recombination pos- 
sibly impacting the rate of evolution of new species. 
Unfortunately, only a single monospore culture of F. 
manipularis was recovered from the field, and other mate- 
rial was not available in this study to determine the nuclear 
status of the mycelium, or to study nonself confrontation 
reactions in mycelial contact zones between different 
dikaryons. In the current study spores from different basi- 
diomata of Filoboletus occurring in proximity could have 
contaminated our samples and the possibility of cloning 
errors and multiple gene copies cannot be completely 
discounted. Most significantly there is evidence from full 
genome sequences for the occurrence of multiple copies of 
the tefla gene in approximately 10% of basidiomycete 
species with no information on the possible occurrence 
of polyploidy or aneuploidy in Filoboletus or close rela- 
tives. Therefore, additional research is necessary to con- 
firm a perhaps common phenomenon of three or more 
component nuclei representing multiple parental monokar- 
yons contributing to the development of individual 
basidiomata. 
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